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Roles of Coulomb interaction, orbital degeneracy and Jahn- Teller coupling in double-exchange 
models are examined for Mn perovskite oxides. We study the undoped Mott insulator as well 
as metal-insulator transitions by hole doping, and especially strong incoherence of ferromagnetic 
metal. We derive models where all the spins are fully polarized in two-dimensional planes as in 
the experimental indications, and investigate their ground-state properties by quantum Monte 
Carlo method. At half filling where the number of e g electron is one per site on average, the 
Coulomb interaction opens a Mott gap and induces a staggered orbital ordering. The opening 
of the Mott gap is, however, substantially slower than the mean-field results if the Jahn- Teller 
coupling is absent. The synergy between the strong correlation and the Jahn- Teller coupling 
largely enhances the Mott gap amplitude and reproduces realistic amplitudes and stabilization 
energy of the Jahn- Teller distortion. Upon doping, the orbital ordering stabilized by the Coulomb 
interaction is destroyed immediately. Toward the metal-insulator transition, the short-ranged 
orbital correlation is critically enhanced in metals, which should be related to strong incoherence 
of charge dynamics observed in experiments. Our model, moreover, exhibits a uniform ordering 
of d x 2_ y i orbital in a wide region of doping in agreement with experimental indications. 

KEYWORDS: perovskite manganites, double exchange model, Coulomb interaction, orbital degeneracy, Jahn- 
Teller coupling, Mott gap, orbital ordering, Jahn- Teller distortion, metal-insulator transition, 
ferromagnetic metal, incoherent charge dymanics, quantum Monte Carlo method 



'Simple' double-exchange (DE) models, that is, mod- 
els with a single non-interacting conduction band fer- 
romagnetically counkjd n wtb localized spins, have been 
intensively studiedBodQEP to understand physical prop- 
erties of Mn pfimvskite oxides, especially collosal mag- 
netoresistanceBElBBEP These models have successfully 
explained several experimental aspects at least qualita- 
tively; for instance, ferromagnetic metals in a hole-doped 
region with a transition to a paramagnetic state by rais- 
ing temperature, and a large negative magnetoresistance 
near the transition. All these properties are due to the 
so-called double-exchange mechanism; doped holes tend 
to gaj^ kinetic energy by aligning localized spins in par- 

However, many open problems still remain. In par- 
ticular, we focus here the following: (i) All the mate- 
rials are the Mott insulator with a charge gap of order 
of eV at half filling (one e g electron per Mn site on av- 
erage), where two-dimensional (2D) anisotropy realizes 
in spin and orbital onjjgdjigs^ pjd_.a cooperative Jahn- 
Teller (JT) distortionMBBB'rJ'Ef (ii) Upon doping of 
holes, a metal-insulator (MI) transition occurs from the 
anisotropic Mott insulator to an isotropic ferroma, 
metal through disordered ferromagnetic insulator.1 
(iii) In the ferromagnetic metallic state, charge dynam- 
ics shows strong incoherence with tiny Drude weight even 
at low temperatures whjlethe linear specific-heat coef- 
ficient 7 remains small.EZHifP All these problems suggest 
importance of elements neglected in the simplified mod- 
els; Coulomb interaction, orbital degrees of freedom, JT 



coupling and so on. 

In several recent studies, effects of these factors previ- 
ously neglected WQ££ also explored; mean-field approxi- 
mation,!^' exact diagonalization of small civ 



ters,E§ Gutzwiller techniqueJ^Lslave-fermion theory,^ 
dynamical mean-field theoryEZrEP and perturbation the- 
ory.Ej' Nevertheless it is still controversial which of the 
electron correlation, orbital degeneracy or the JT cou- 
pling is the driving mechanism of experimental proper- 
ties, particularly strongly incoherent charge dynamics, 
and what are realistic parameters of theoretical models. 
These need to fully consider quantum fluctuations. 

In this work, we investigate DE models including 
Coulomb interaction, orbital degeneracy of conduction 
bands and JT coupling. Especially, we discuss quantum 
fluctuation effects of on-site Coulomb interaction on half- 
filled states as well as on MI transitions by hole doping 
in the ground state. In order to obtain unbiased results 
under strong quantum fluctuation, we. .use here quan- 
tum Monte Carlo (QMC) techniqueMeJ Our results 
at half filling show that a strong Coulomb interaction 
is required to understand undoped manganites au 
tatively. Experimentally-observed JT distortionaZrEl 
are reproduced as a consequence of the synergy between 
the strong correlation and the JT coupling. Nature of the 
MI transition is investigated in detail; we show that this 
transition is a continuous one with a critical enhance- 
ment of orbital correlations—wldch may be relevant to 
incoherent charge dynamics. EZl'tfP A uniform orbital or- 
dering observed in doped materialsE 2 ]' Ei* is also repro- 
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duced in our model. 

The model we study here consists of three terms as 



H — Hd + W c i_ph + H p h 



(1) 



The first term is derived in a particular limit of DE 
models with on-site Coulomb interaction as well as with 
twofold degeneracy of e g orbitals. The original form 
of Ti c \ is given by Hbe = H t + Hint + Hk, where 
n t = T,ij,vv>,<, t ij'4v<T C j«'<* is the hopping term, here 
Civo{c\ U(7 ) annihilates (creates) a cx-spin e g electron at 
site i = 1, • • •, Ns with orbital v = 1 or 2; H.- m t contains 
on-site interactions within the doubly-degenerate e g or- 
bitals; and Hr — K ^i" ■ S* 2 " denotes the Hund's rule 
coupling between e g electmns. and t% g localized spins. 
Following previous studiesp'tm'tP we consider the limit 
of A'|S t2 «| > t\f . In this limit, aligns parallel to 
S/ 9 in each site, therefore is renormalized by a 
relative angle of ti g spins in sites i and j. Moreover, 
we assume perfect polarization of spins in the ground 
state based on the following arguments. At finite hole 
doping, it is known that 7Yde shows a metallic ground 
state with perfect polarization of spin^ whirh rnrre- 
sponds well to many doped manganites.^EZfEa'Eil'H At 
half filling, these materials commonly become an A-type 
antiferromagnetic insulator which consists^oLfexromag- 
netic layers stacking antiferromagnetically.lDOEif These 
show that the 2D structure of perfectly-polarized ferro- 
magnetic planes may persist in both sides of MI tran- 
sition caused by hole doping. Therefore, within a 2D 
plane at low temperature, spin degrees of freedom may 
be dropped. Then the Hamiltonian in the plane reads 



Wei = £ i W 4< 
ij,vu' 



5X^-5), (2) 



where rij„ = c iv Ci V is a number operator. The transfer 
matrices are renormalized independently of sites because 
of the perfect polarization; and the interaction U is not 
a bare interorbital Coulomb repulsion U12 but U12 — J12, 
where J12 is the Hund's rule coupling between orbitals 
v = 1 and 2. Here hopping integrals between nearest- 
neighbor sites in the 2D plane are explicitly given as 
t 11 = -3t /4, t 22 - -to/4, t 12 = t 21 = -(+)v% /4 in 



the x(y) r 



irection for d x 2_ y i{y = 1) and d 3z 2 _ r i(y = 2) 



orbitals,^ where to is defined later. This model can be 
viewed as a generalization of 'ordinary' Hubbard models 
where t vu oc 5 vu i In the following, we consider the 
above 2D model partly because of its numerical feasibil- 
ity. Although it ignores some of important effects due 
to three-dimensionality, it may contain at least nontriv- 
ial important physics of the MI transition by hole doping 
where the 2D anisotropy was suggested to be relevant .Ell* 
Moreover, this model is useful to discuss an orbital or- 
dering at half filling since it also has 2D structure as the 
spin ordering. We also note that our model may have 
some relevance to layered Mn compounds,E*P which con- 
sists of stacking of single or double Mn02 layers with 
rather small crystal-field splitting of d x 2_ y 2 and d 3z 2_ r 2. 

The rest terms in eq_M) denote the JT coupling 
in a classical treatment .ESIHj) These should be appro- 
priate when lattice distortion has a long-ranged order. 



In perovskite manganites, a cooperative JT distortion 
which also has 2D anisotr-GpY.-as spin and orbital order- 
ing appears at half filling. Qllil'lla' Because of the symme- 
try of e- .w^jve functions, we consider only two phonon 
modes.cj'liyE 2 ]) Then we obtain the forms as 



(3) 



where g and I are two-component vectors defined as 
Si = -2^(008 20;, sin 20*) and I 4 = (3?,3f); g and 
k are the electron-phonon coupling and the spring con- 
stant. Here Ui denotes displacement of oxygens sur- 
rounding the i-th Mn site and 9i determines the coupling 
angle. T]f*(/x = x,y,z) is called the pseudo-spin opera- 
tor which denotes two degrees of freedom of orbitals as 
= \ ~YLw' ^vv'4v c iv' ^ where a is the Pauli matrix. 
When the oxygen between the i and j-th Mn sites 
moves by Ui, the hopping integral between Mn sites 
changes since originally it consists of a product of over- 
lap integrals between Mn and O sites. Note that in both 
phonon modes, oxygens shift only along the direction 
of Mn-O-Mn bonds. A displacement of Ui .modifies the 
Mn-0 hoppings by a factor of (1 ± Ui) 7 ^ 2 ^ therefore 
we define a unit of the Mn-Mn hopping integrals in eq. 

d) as < =to(l + ^) 7/2 (l-^) 7/2 - 

In the following, we discuss ground state properties of 
the model eq. _(3Lin two dimensions by the projection 
QMC method.tSEIP Computational details including the 
negative sign problem will be described elsewhere.E3> We 
take energy and length units as to and the lattice con- 
stant, respectively. From experiments and band calcu- 
lations, to is considered ijQpiiKL^muAd-p-S °V and U is 
estimated at several eV.EjEJIlJ'EJ'OEZP In the present 
study, we change U /to as a parameter and discuss its 
realistic value based on our results. The spring constant 
k is roughly estimated by the frequency of anpOxygen 
bond stretching phonon at order of 10 or 100 eVO The 
electron-phonon coupling g is not easy to determine ex- 
perimentally, however, it may be several eVO We set 
here k = 100 and g = 10. The parameters Ui and Qi are 
determined by a mean-field scheme for simplicity; we fix 
9i at experimental values as mentioned below and deter- 
mine a uniform solution Ui = u to minimize total energy 
obtained from the QMC calculations. 

First, we discuss results at half filling (Xhj, n iv — As). 
Fig. 1 shows the U dependence of charge gap and stag- 
gered moment of orbital polarization in the absence of 
the JT coupling. For the comparison, QMC re 
for 'ordinary' 2D Hubbard models are also plottec 
with the mean-field results. Here, the charge gap is 
calculated from shifts of chemical potential upon dop- 
ingJ-S' In our models, the charge gap opens consider- 
ably slower than the case of 'ordinary' Hubbard models, 
which leads to a remarkable deviation from the mean- 
field results, presumably because of their characteris- 
tic nesting propertycP In order to reproduce the real- 
istic charge-gap amplitude within the electronic model 
eq. d), we need a much stronger interaction. Fig. 1 
(b) shows that the interaction U induces the staggered 
long-ranged orbital ordering of T x component, that is, 
(\d x 2_ y 2) + \d 3z 2_ r 2))/(\d x 2_ y 2) - \d 3z 2_ r 2)) type, whose 
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pattern is slightly different from the experimentally- 
observed one, \d^ x 2_ r 2) /\d^ y 2_ r 2) . This discrepancy has 
also been pointed oiU_.ip_.previous studies by the mean- 
field approximations]^' E3 The moment grows rapidly as 
the values of U, probably as a consequence of its uniaxial 
anisotropy. Our results show that the Coulomb interac- 
tion U by itself leads to a Mott insulating state with a 
staggered orbital ordering, although its ordering pattern 
does not show complete agreement with experiments. 



1 2 3 4 5 




u 

Fig. 1. QMC results (circles) without JT coupling for (a) the 
charge gap and (b) the staggered orbital polarization. The data 
are shown in the limit of N§ — > oo, which are obtained by size 
extrapolations from jVg = 4 X 4 to 10r*i Wi Squares are QMC 
results for 'ordinary' Hubbard modelscj'EJ The gray (dotted) 
curve shows the mean-field results for our model (for the 'ordi- 
nary' Hubbard model). 

If we take the JT coupling into account, importance 
of electron correlation becomes clearer. Fig. 2 shows 
ground state energy per site as function of oxygen dis- 
placement u. Here we assume 9i = (— l)l r r+ r « >n/6, fol- 
lowing the JT pattern observed in experiments. The 
total energy is substantially lowered and minimized at a 
finite u; the stabilization energy Ejt — 0.13 at u ~ 0.04 
for U = 4. At the same time, the orbital-ordering struc- 
ture changes to the experimentally-observed one, namely 
|d 32 ,2_ r 2)/|d3j ) 2_,.2). In addition, the charge gap ampli- 
tudes become considerably larger than those without the 
JT coupling; for instance, at U = 4, A c 1 (~ 1.1 for 
hole doping and 1.7 for electron doping). These stabi- 
lization energy and charge gap grow with increasing U 
or g/k. This indicates suppression of quantum fluctu- 
ations as a cooperative and combined effect of U and 
g/k. Experimentally, u is about 0.04 and Ejt is below 
or around 0.1 eV in the light of the structural transition 
temperature (for instance, ~ 800K for LaMn0 3 )MEl>B> 
It is quite surprising that in spite of neglecting the 
three-dimensionality, our present results show quantita- 
tive agreement with these experimental results; although 
more tuning of U or g/k (probably larger U) is required 
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to obtain more realistic value of the Mott gap. We should 
stress that in the absence of U, the JT distortion is far in- 
sufficient in reproducing these results. Synergetic effects 
of the strong correlation and the JT coupling play a cru- 
cial role to reproduce the realistic JT distortion, orbital 
order configurations and the Mott gap amplitude. H 
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Fig. 2. Ground state energy per site at half filling as function of 
the oxygen distortion u. We take U = 4, g = 10 and k = 100. 
The curve is a polynomial fit to the data at TVg = oo. 
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Fig. 3. Inverse plot of peak values of the orbital correlation func- 
tion as function of doping concentration at U = 4. The gray line 
is the linear fit to the data at < <5 < 0.15. 

Doped holes lead to an MI transition from this Mott 
insulating state. We focus effects of Coulomb interac- 
tion on this transition. The ferromagnetic insulator pre- 
sumably with the JT distortion adjacent to the Mott 
insulator is beyond the scope of this paper because the 
JT coupling was not studied in detail in doped cases. 
Our data of doping concentration dependence of chemi- 
cal potential indicate that there is no phase separation. 
In Fig. 3, we plot the inverse of peak values of the 
orbital correlation function T X (Q) where Q = {%, ii) 
as function of hole concentration 5. Here we define 
T"(k) = Wg^EyZfT/expCikry) (/i - x,y,z). Our 
results show that the staggered orbital correlations ap- 
pear to be insensitive to the system sizes for 5 > 0. More- 
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over, for small S, they are scaled as T^Q) -1 ~ S. These 
indicate that the staggered orbital polarization stabilized 
by U is destroyed simultaneously with the MI transition 
at 5 = 0. In 'ordinary' Hubbard models, the similar scal- 
ing of spin correlations and an anomalous scaling of com- 
pressibility were found.© These suggest a novel univer- 
sality class of the MI transition where residual entropy of 
spin degrees of freedom may induce, strong incoherence 
of a metal near the critical point E3 In our model, al- 
though we have not explicitly studied the compressibility 
anomaly yet, the critical behavior of orbital correlations 
and its implication of large entropy may well be tlp.oxigin 
of a strong incoherence observed in experiments. 

In doped cases_-Qur model exhibits a uniform ordering 
of d x 2_ y 2 nrhit.a.l N) which has been suggested in recent 
experiments.Ej'Ej Our results show that a positive peak 
of T z (k = 0) grows as S and has a maximum around 
S = 0.5. This ordering presumably comes from an op- 
timization of the kinetic energy in the 2D planes, since 
it is observed even in the non-interacting 2D model. We 
note that this ordering is substantially enhanced by U . 

To summarize, we have investigated ground states of 
double-exchange models with on-site Coulomb interac- 
tion, twofold degeneracy of conduction bands and JT 
coupling using the QMC method under the condition of 
perfect-polarization of spins in 2D planes. Experimental 
results of undoped manganites, especially the JT distor- 
tion, are quantitatively reproduced with realistic param- 
eters. Importance of the synergy between the strong elec- 
tron correlation and the JT coupling is clarified. Upon 
doping, our models show a continuous metal-insulator 
transition with a critical enhancement of staggered or- 
bital correlations. We have speculated its relation to 
strong incoherence of charge dynamics in ferromagnetic 
metals observed in experiments. In addition, our mod- 
els exhibit a uniform ordering of d x 2_ y 2 orbital in doped 
cases, which agrees with experimental indications. 
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